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Abstract A hydrodynamic approach is used to calculate an asymptotics of the Emptiness 
Formation Probability - the probability of a formation of an empty space in 
the ground state of a quantum one-dimensional many body system. Quantum 
hydrodynamics of a system is represented as a Euclidian path integral over con- 
figurations of hydrodynamic variables. In the limit of a large size of the empty 
space, the probability is dominated by an instanton configuration, and the prob- 
lem is reduced to the finding of an instanton solution of classical hydrodynamic 
equations. After establishing a general formalism, we carry out this calcula- 
tion for several simple systems - free fermions with an arbitrary dispersion and 
Calogero-Sutherland model. For these systems we confirm the obtained results 
by comparison with exact results known in Random Matrix theory. We argue 
that the nonlinear hydrodynamic approach might be useful even in cases where 
the linearized hydrodynamics fails. 
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1. Introduction 

The hydrodynamic approach to correlation functions in quantum many body 
systems has a long history [1, 2]. Generally, hydrodynamic equations are non- 
linear and dispersive. However, usually, it is the linearized hydrodynamics 
which is used to extract long distance asymptotics of correlation functions [2]. 
In particular, in one spatial dimension the linearized quantum hydrodynamics 
(or bosonization) is especially useful [3]. It was argued that corrections to the 
bosonization due to the nonlinearities of the hydrodynamics (curvature of an 
underlying fermionic spectrum) are irrelevant in calculation of leading terms 
of correlation functions [4]. 

Although the linearized hydrodynamics turned out to be a very powerful 
tool in studies of correlation functions, there are many phenomena for which 
nonlinearities are essential. The goal of this lecture is to consider a particular 
example where the linear approximation fails but nonlinear hydrodynamics can 
still be used to extract non-perturbative results for correlation functions. 

As such an example we consider a particular quantity which plays an impor- 
tant role in the theory of quantum, one-dimensional, integrable systems [5] - 
the Emptiness Formation Probability (EFP). This correlation function was ar- 
gued to be the simplest of correlators in some integrable models [5]. The EFP 
P(R) is essentially, a probability of formation of an empty region of a size 2R 
in the ground state of the many body system. In integrable models the EFP has 
an exact representation in terms of determinants of Fredholm operators [5, 6] 
or multiple integrals [7]. These expressions are exact but are very complex, 
and extracting, e.g., long distance asymptotics R — ► oo from the exact expres- 
sions is a non-trivial problem (see Appendix A and references therein for some 
known results). However, the limit of large R is precisely the limit where hy- 
drodynamic description is applicable, and we illustrate that it, indeed, correctly 
produces the leading term of the asymptotics of the EFP. The hydrodynamic 
approach can also be used to calculate more complicated quantities than EFP 
(which is, essentially, not a dynamic but a ground state property). Moreover, 
the applicability of hydrodynamics is not limited to integrable systems. 

We start with some definitions as well as with making a connection to Ran- 
dom Matrices. 

Emptiness Formation Probability. Consider a one-dimensional quantum 
liquid at zero temperature. The wave function of the ground state of the liquid 
*g(xi, X2, ... ,xn) gives the probability distribution \^g\ 2 of having all N 
particles at given positions xj, where j = 1, . . . , N. We introduce the Empti- 
ness Formation Probability (EFP) P(R) as a probability of having no particles 
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with coordinates —R<xj<R. Formally we define 

P(R) = L v / dx 1 ...dx N \V G (x 1 ,...,x N )\ 2 , (1.1) 

\*G|^G/ J\xj\>R 

or following Ref . [5] 

P(i?) = lim (^ G |e- Q ^-V (:E)cte |^ G ), (1.2) 

where p{x) is a particle density operator 

N 

p( x ) = ^2 s ( x - x j)- ( L3 ) 

We are interested in an asymptotic behavior of P(R) as pqR — > oo, where 
po is the average density of particles in the system. EFP -P(-R) gives us the 
probability that the one-dimensional river parts to make a ford of a macroscopic 
size 2R. 

Random Matrices. The EFP (1.1) introduced for a general one-dimensional 
quantum liquid is a well-known quantity in the context of spectra of random 
matrices [8]. Namely, it is the probability of having no eigenvalues in some 
range. Consider e.g., the joint eigenvalue distribution for the Circular Unitary 
Ensemble (CUE). The CUE is defined as an ensemble of N x N unitary ma- 
trices with the measure given by de Haar measure. Diagonalizing matrices and 
integrating out unitary rotations, one obtains [9] 



jDu^jUde, n 



\ e i9j _ e i9 k 



j=l l<j<k<N 



(1.4) 



where (3 = 2 for CUE and e i9 * with j = 1,...,N are the eigenvalues of a 
unitary matrix. One can read the joint eigenvalue distribution 



P N (9i, ...,9n) = const. Jl \e ie i - e 

l<j<k<N 



lti k 



(1.5) 



Here we replaced the sums over particles (eigenvalues) to integrals with par- 
ticle densities. We left only terms which are dominant in the limit of large 
N. 1 Now we introduce the probability of having no eigenvalues on the arc 

'In fact, one can do a better job including subdominant corrections. See Ref.[9] for details. 
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—a < 9 < a as 



" J erf [-«,<*] fj[ 




(1.6) 



This quantity is known as ^(O, a) in notations of [9]. For orthogonal, unitary, 
and symplectic circular ensembles the joint eigenvalue distribution is given by 
(1.5) with [3 = 1, 2, 4 respectively. 

Spin chains and lattice fermions. The EFP can also be defined for spin 
chains where we are interested in the probability of having a ferromagnetic 
string of the length n in the ground state of the spin chain. The Jordan- Wigner 
transformation maps spin 1/2 chain to a one-dimensional lattice gas of spinless 
fermions. Under this mapping the ferromagnetic string corresponds to a string 
of empty lattice sites and one can write 



where ipj , ijr- are annihilation and creation operators of spinless fermions on 
the lattice site j. Therefore, the probability of the formation of a ferromag- 
netic string in spin chains corresponds to the Emptiness Formation Probability 
of Jordan-Wigner fermions. We are going to use a language of particles in 
this paper but all results are also valid for corresponding one-dimensional spin 
systems. 

2. Instanton or rare fluctuation method 

In the limit of large R (poR 3> 1) we use a collective description instead 
of dealing with individual particles. We assume that u is some collective field 
describing the hydrodynamic motion of a one-dimensional liquid. 2 The dy- 
namics of the liquid is defined by a Euclidian partition function 




(1.7) 




(2.1) 



where S[u] is the Euclidian action 




1/2T 



(2.2) 



2 Later we will use the conventional displacement field as u. See Eq. (3.3). 
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and the inverse temperature 1/T defines periodic boundary conditions in the 
imaginary time r. 

The asymptotic behavior of P(R) in the limit of large R is defined by a rare 
fluctuation when all particles move away so that at some time t = we have 
no particles in the spatial interval [-R, R]. Then with an exponential accuracy 

P(R) ~ e~ Sopt . (2.3) 

Here S op t is the value of the action (2.2) on the trajectory u(x, t) which min- 
imizes (2.2) and is subject to EFP boundary conditions. These are: the EFP 
boundary condition 

pit = 0; -R < x < R) = 0, (2.4) 
and standard boundary conditions at infinity 

P Po, x,r^oo, 

v — > 0, x,r^oo, (2.5) 

where p, u are the density and velocity of particles related somehow to the 
collective field u. 

Let us estimate the probability of such rare fluctuation 3 . Without the bound- 
ary condition (2.4) the minimum of the action subject to (2.5) is obviously 
given by the constant solution p(x, t) = po, v(x, t) = 0. The condition (2.4), 
however, disturbs this constant solution in some area of space-time around the 
origin. The spatial extent of this disturbance is of the order of R. If we assume 
that the quantum system we are dealing with is some compressible liquid we 
expect that the typical temporal scale of the disturbance is of the order of R/v s , 
where v s is a sound velocity at p = po. We conclude that the (space-time) 
"area" of the disturbance scales as R 2 and the action S op t ~ R 2 . Therefore, 
we expect a Gaussian decay for (2.3) 

P(R)~e~ aR2 , (2.6) 

where a is some (non-universal) constant depending on the details of (2.2). 

This argument can be extended for the case of low but finite temperature. 
Namely, while the temporal extent of the instanton R/v s is smaller than the 
inverse temperature R/v s <C 1/T instanton "does not know" that the temper- 
ature is not zero and one obtains an intermediate Gaussian decay of EFP (2.6). 
However, for long enough R it is the 1/T scale that defines the temporal size 
of the instanton, the space-time area of the disturbance scales as R and one 
obtains 

P(R) ~ e~ lR . (2.7) 

3 In this section we follow closely the qualitative argument of Ref.[10] 
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One expects a crossover between the Gaussian and the exponential behavior 
taking place at R ~ v s /T. The result (2.7) is very familiar from statistical 
physics with 7 measuring the difference (per unit length) of free energies be- 
tween the true ground state of the liquid and the empty state. The Gaussian 
decay (2.6) is a manifestation of an effective increase of the dimensionality 
from one to two in quantum systems at T = 0. 

In the above argument we used the assumption of compressibility of the 
quantum liquid. For incompressible liquids one has a finite correlation length 
of density fluctuations. The argument can be extended then to obtain an ex- 
ponential form (2.7) with the correlation length playing the role of an inverse 
temperature [11, 12]. 

3. Hydrodynamic approach 

The collective description we are looking for is nothing else but a hydrody- 
namic description in terms of density p and current j = pv (or velocity v ) of a 
one-dimensional liquid. For simplicity, let us consider first the case of a system 
with Galilean invariance. Then we write the Euclidian action of a liquid as 



S = dfx 



f p + pe(p) + ... 



(3.1) 



where the first term is fixed by the Galilean invariance and is the kinetic energy 
of the liquid moving as a whole. The second term is the internal energy of the 
fluid which is determined by the equation of state of the liquid. e(p) is the 
internal energy per particle at given density p. The terms denoted by dots are 
the terms which depend on density and its spatial derivatives. These terms 
will be small in the problems, where density and velocity gradients are small 
compared to p . Let us now remember that due to the particle conservation the 
density and current are not independent variables but are related by a constraint 
- the continuity equation 

d t p + d x j = 0. (3.2) 

One can easily solve the constraint (3.2) introducing a particle displacement 
field u such that 

p = po + d x u, 

j = -d t u. (3.3) 

Microscopic definition of the displacement field is u(x)+pox = Y^j 0(x—Xj), 
where 0(x) is a step function and Xj-s are coordinates of particles. It is easy 
to check that the configuration u(x,t) minimizing S[u] from (3.1) is given by 
5S = or after simple algebra 

d t v + vd x v = d x d p [pe(p)], (3.4) 
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which is the Euler equation of a one-dimensional hydrodynamics. The sign of 
the r.h.s. of (3.4) differs from the conventional minus sign because we work in 
the Euclidian formulation. 

The action (3.1) with the parameterization (3.3) provide us with a varia- 
tional formulation of one-dimensional classical hydrodynamics. To calculate 
the probability of hydrodynamic fluctuations at zero temperature we have to 
quantize this hydrodynamics. To the best of my knowledge, the first "quanti- 
zation" of hydrodynamics was done by L. D. Landau in Ref.[l], where he used 
essentially (3. 1) as a quantum Hamiltonian of the liquid with density and veloc- 
ity fields satisfying commutation relations [p(x),v(y)] = —id x 5(x — y). For 
the purpose of evaluating a rare fluctuation the path integral approach is more 
useful and we use the partition function (2.1), where u(x,t) is the displace- 
ment field and the action S[u] is given by (3.1) and (3.3). We notice here that 
we did not specify the measure of integration Du in (2.1). Finding this mea- 
sure requires a derivation of an effective hydrodynamic formulation from the 
underlying microscopic physics. However, the "non-flatness" of the measure 
gives only gradient corrections similar to the already omitted terms denoted 
by dots in (3.1). Those corrections to the measure and the action will produce 
subdominant contributions to the asymptotics of EFP and will be neglected in 
this work. 

Let us now summarize our strategy for the calculation of the leading term 
of the EFP. We solve classical equations of motion (3.2,3.4) with EFP bound- 
ary conditions (2.4,2.5) and then find the value S op t of (3.1) on the obtained 
solution. Finally, (2.3) will give us the dominant contribution to the EFP at 

R — > oo. 

We conclude this section with two remarks. First, it will be convenient for 
us to generalize the problem and replace the EFP boundary condition (2.4) 
by a slightly more general depletion formation probability (DFP) boundary 
condition [10] 

p(t = 0;-R<x<R) = p, (3.5) 

where p is some constant density. In the case p = we obtain the EFP problem 
while for p close to po one can use the bosonization technique to calculate 

P(R;p). 

Second remark is that one can easily obtain the functional dependence of 
P(R) even without solving hydrodynamic equations (3.2,3.4) using very sim- 
ple scaling arguments. Indeed, the equations (3.2,3.4) are uniform in space 
and time derivatives so that if p(x, t) and v(x, t) are solutions, then p(Xx, At), 
v(Xx, At) are also solutions of hydrodynamic equations. Choosing A = R we 
obtain the boundary condition (3.5) as p(t = 0; — 1 < x < 1) = p and the 
only dependence on R is left over in the integration measure of (3.1), which 
gives S op t ~ R 2 - Thus, we obtain a Gaussian decay of DFP (or EFP in partic- 
ular) as a function of R. Corrections to this behavior come from the terms of 
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higher order in gradients in the hydrodynamic action (denoted by dots in (3.1)) 
as well as from fluctuations around the saddle point (classical) trajectory in the 
partition function (2.1). 

4. Linearized hydrodynamics or bosonization 

Before proceeding to the general case, let us consider the DFP for 

^ « 1, (4.1) 

Po 

i.e., the probability of formation of a small constant density depletion along 
the long string — R < x < R. In this case the deviation of the density from po 
is small almost everywhere (see below) and one can use a linearized version of 
generally nonlinear hydrodynamic equations (3.2,3.4). Expanding the classical 
action (3.1) in gradients of the displacement field u (3.3) we obtain in harmonic 
approximation 4 

S bos = — [ d?x\{d,uf, (4.2) 
Po J * 

where we scaled v s ot — > t. The sound velocity v s q at p = po is defined as 



v % = P d l(P e (p)) 



(4.3) 



P=P0 

The corresponding equation of motion is the Laplace equation 

Au = (4.4) 

with the DFP boundary condition 

u(x, t = 0) = -(p - P)x, for - R < x < R. (4.5) 

It is easy to see that the solution of (4.4,4.5) decaying sufficiently fast at 
infinity is given by 



u(x,t) = -(p - p) Re \ zo -^zl-R?y (4.6) 

where we introduced the complex notation zq = x + iv s ot (where t is the orig- 
inal imaginary time). Indeed, at i = 0, — R < x < R the complex coordinate 
zo is real and square root in (4.6) is purely imaginary so that (4.5) is satisfied. 
At space-time infinity zq — > oo we have 

u(x,t) « -— - (4.7) 
zo Zo 



4 The harmonic approximation to the nonlinear hydrodynamics of quantum liquid is equivalent to a linear 
bosonization approach to interacting one-dimensional particles. 



10 



with 



a = a 



= -^-p)R\ 



(4.8) 



We obtain from (4.7) that at zq — ► oo 



a a 



P ~ P0 H 2+^2' 



•o ^0 




(4.9) 



which obviously satisfy the boundary conditions at infinity (2.5). 

Now that we obtained the solution of hydrodynamic equations it is a straight- 
forward problem to calculate the value of the action (4.2) on this solution 5 and 
obtain 



We note that the linearization of hydrodynamic action (4.2) is not self- 
consistent near the ends of the string t = 0, x = ±R where the solution 
(4.6) is singular and gradients of u diverge. However, one can estimate the 
corrections coming from those areas and find that they change the coefficient 
[tt(Po — p)\ m front of the R 2 in (4.10) by terms of the higher order in small 
parameter (4.1) [10]. 

We also note here that although the value of the coefficient a given by (4.8) 
is approximate and is valid only in the limit of a very weak depletion, the 
asymptotic forms (4.7,4.9) are very general as they depend only on the lin- 
earization of hydrodynamic equations at x, t — > oo, where it is always possi- 



5. EFP through an asymptotics of the solution 

The calculation of the EFP has already been reduced to the calculation of the 
value of the classical action on the solution of equations of motion satisfying 
EFP boundary conditions. In this section we are going use a Maupertui prin- 
ciple [13] to obtain a simple expression for S op t in terms of the asymptotics of 
the EFP solution of hydrodynamic equations. 

Let us calculate the variation of the action (3.1) with respect to the displace- 
ment field u 



S DF p = ^Hpo-P)R?- 

Z IT On 



(4.10) 



ble. 



5S 




(5.1) 



5 See the next section on how to avoid doing this calculation. 
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We kept here surface terms (full derivatives) in addition to the last term which 
produces the equation of motion (3.4). Now we assume that the action S is 
calculated on the EFP (or DFP) solution of equations of motion and consider 
the derivative of this action with respect to the equilibrium background density 
Pq. We have 



-d t (vd Po u) - d x 



u 



(5.2) 



where we used the fact that u satisfies the Euler equation and dropped the 
last term in (5.1). (5.2) has only derivative terms and can be re-written as a 
boundary contribution 



j> vd Po udx + (d p (ep) - d po udt 



(5.3) 



where the integral is taken over the infinitely large contour around xt plane. 
As the integrand in (5.3) should be calculated at (infinitely) large x and t we 
can use in (5.3) the general asymptotics (4.7,4.9). After simple manipulations 
we obtain our main result 



d P0 S O pt = 2tt— (a + a). 
Po 



(5.4) 



As a simple check of this result we take the value of a obtained in bosoniza- 
tion approach (4.8) and substitute it into (5.4). We immediately obtain d Po S op t = 
^—(po — P)R 2 , which is equivalent to (4.10) up to the terms of higher order 



PQ 



in the small parameter (4.1). 



6. 



Free fermions 



In this section we find the EFP for a free Fermi gas in one dimension. First, 
let us find the internal energy e(p) of the gas. The density of fermions is given 
in terms of Fermi momentum kp as 



kp 



dk 

2tt 



kp 

7T 



The energy per unit length is 

pAp) = 



1,3 
ivp 



7T- 



(6.1) 



(6.2) 



where we put h = 1 and fermion mass m = 1. The energy per particle in a 
free Fermi gas with density p is 



/x 71-2 2 
<P) = -JP- 



(6.3) 
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We calculate the sound velocity using 

v 2 s = P d 2 p {pe) (6.4) 

and obtain 

v s = irp = k,F (6.5) 

- the well known result that the sound velocity and Fermi velocity are the same 
(remember that in our notations m = 1 and vp = fcp)- 

Hydrodynamic equations (3.2) and (3.4) for free fermions are 

d t p + d x (pv) = 0, 
d t v + vd x v = n 2 pd x p. (6.6) 

Introducing a complex field 

w = up + iv, (6.7) 

we re-write both equations (6.6) as a single complex Hopf equation 6 

dtw — iwd x w = 0. (6.8) 

The latter equation has a general solution 

w = F(z), (6.9) 

where F(z) is an arbitrary analytic function of a complex variable z defined as 

z = x + iwt. (6.10) 

The boundary conditions will determine a particular analytic function F(z) so 
that the equations (6.9,6.10) will define the solution w(x, t) of (6.8) implicitly. 
It is easy to check that the unknown function F(z) for DFP is given by 

F{z) =irp + n(p -p)-j=L= (6.11) 

and the one for EFP can be obtained by putting p = in (6.1 1). The space-time 
configuration of the density and velocity fields minimizing the hydrodynamic 
action is given by real and imaginary parts of w(x, t) (see Eq. (6.7)) which is 
implicitly defined by 

w-Tip = +7r(/9 - p) ^ 2 _ R2 (6- 12) 

with (6.10). This solution is relatively complicated (see Figs. 1, 2) and a direct 



6 In real time formalism, instead of w, w one introduces "right and left Fermi momenta" kn ^ = 7rp ± v 
which satisfy the Euler-Hopf equations d t k + kd x k = reflecting the absence of interactions between 
fermions. 
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Figure 1. The density profile p(x, t) is shown for the EFP instanton. The density diverges at 
points (x, t) = (±ii, 0). The shape of the "Emptiness" is shown in Fig.2. 



calculation of the value of the hydrodynamic action (3.1) on this solution is 
cumbersome. However, one can get the result very quickly using (5.4). Indeed, 
in the limit x, t — > oo we have w — ► irpo and z = x + iwt — » x + iirpot = Zq. 
Therefore, the asymptotics of (6.12) is given by 

R 2 

w - vrpo ~ vr(p - p)tT2 ( 6 - 13 ) 
Zz 

and taking, e.g., its real part 

P ~ Po ~ \(Po ~ P)R 2 ( 4 + 4) • (6-14) 
Comparing (6.14) with (4.9) we extract 

a = ^(/'o-/o) J R 2 (6-15) 
and substituting into (5.4) we derive d Po S op t = ir 2 (po — p)R 2 and 

Sopt = ~ k(po - P)^] 2 • (6-16) 
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Figure 2. The region of the x — t plane in which p(x, t) = for the EFP instanton for free 
fermions is shown. The boundary of the region is given by an astroid x' 2 ^ 3 + (irpot) 2 / 3 = R 2 ^ 3 . 

This gives for DFP and EFP probabilities respectively 

Pdfp(R) ~ exp|-i[^(p -p)i?] 2 |, (6.17) 

Pefp(R) ~ exp j-^vrpoi?) 2 }. (6.18) 

In the case of free fermions an exact asymptotic expansion of EFP in l/R 
is known [8] with first few terms given by (B.2). The instanton contribution 
(6.18) gives an exact first (Gaussian) term of this expansion. 

7. Calogero-Sutherland model 

Our next example of a one-dimensional liquid is the Calogero-Sutherland 
model - a model of one-dimensional particles interacting with an inverse square 
potential. The Hamiltonian of the model is 



k=l,k^i 

Here we again use the units h = 1 and m = 1. This model is known to be 
integrable [14, 15]. We are interested in the thermodynamic limit of (7.1). The 
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easiest way to go to this limit is to consider (7.1) in an additional harmonic 
potential which does not destroy the integrability of the model and then take a 
limit of number of particles N —> oo and the strength of the potential going to 
zero so that the density is kept constant and equal po. We omit all these details 
which can be found in the original papers [14, 15]. We mention only that the 
ground state wave function of (7.1) is 

^GS = H(x j -x k ) x (7.2) 

j<k 

and shows an intermediate statistics interpolating between non-interacting bosons 
(A = 0) and non-interacting fermions (A = 1). We also notice here that the 
probability distribution of particle coordinates 

\*gs\ 2 = U\x,-x k \ 2X (7.3) 

j<k 

at particular values of coupling constant A = 1/2, 1, 2 coincides with the joint 
probability of eigenvalues for the orthogonal, unitary, and symplectic random 
matrix ensembles respectively (see Eq. (1.5)). 

To calculate the leading behavior of the EFP for the Calogero-Sutherland 
model we need the internal energy e(p) which can be easily found [15] 

e(p) = yAV . (7.4) 

The (7.4) differs from the free fermion case (6.3) by a factor of A 2 and coin- 
cides with the latter (as expected) at A = 1. Introducing 

w = Xirp + iv (7.5) 

and repeating the calculations of the previous section we obtain 

Pdfp(R) ~ exp^-h[7r(p -p)R] 2 Y (7.6) 

Pefp(R) ~ exp|-iA(7r / 9 0j R) 2 |. (7.7) 

Comparing to the known exact result (B.4) we see that (7.7), indeed, gives 
the exact leading asymptotics of the EFP for the Calogero-Sutherland model. 
Subleading (in l/R) corrections to (7.7) are due to gradient corrections to the 
hydrodynamic action (7.4) and to quantum fluctuations around the found in- 
stanton. 



16 



8. Free fermions on the lattice 

The goal of this section is to illustrate that the hydrodynamic method we 
used is not limited to Galilean invariant systems. We use the method to calcu- 
late the EFP for a system of non-interacting lattice fermions with an arbitrary 
dispersion e{k). Semiclassically, one can describe the evolution of the degen- 
erate ID Fermi gas in terms of two smooth and slow functions k R ^(x,t) - 
right and left Fermi points respectively. The equations of motion are given by 

d t k R , L + e'{k R:L )d x k R , L = 0, (8.1) 

which is an obvious consequence of the absence of interactions. Indeed, the 
derivative e'{k R ^L) is the right (or left) Fermi velocity of particles, and (8.1) is 
nothing else but the statement that the momentum of a particle does not change 
in time. The classical action of ID Fermi gas should reproduce (8.1) as well as 
give the correct energy of the system E = f dx f^/^) ^e{k). 

We assume here that the dispersion of fermions is parity invariant e(k) = 
e(—k) and consider the following Euclidian action 

S = J d 2 x j^- [wd^dtw - wd^dtw] + J • (8-2) 

Here e(k) is the dispersion (including chemical potential) of free fermions 
and w, w are independent variables which are Euclidian versions of k R , — 
We identify the density of particles as p = The equations of motion 

corresponding to (8.2) 

idtw + d x e{w) = 0, 
-id t w + d x e{w) = (8.3) 

are, indeed, the Eucl idian versions of (8.1). Equations (8.3) are complex conju- 
gates of each other (e(w) = e(w) as e(k) is a real function). The difference of 
these two equations give the continuity equation (3.2) with the current j = 
The path integral with the action (8.2) is the integral over two fields w, w. We 
integrate out their difference w — w in a saddle point approximation. Namely, 
we re-write (8.2) in terms of combinations w — w and w + w. Then, we take 
a variation of the action with respect to w — w. This produces the continuity 
equation relating w and w. We solve the continuity equation introducing the 
displacement field. Thus, we have w and w in terms of u as the saddle point 
trajectory. Namely, 

p + d x u, 

-d t u. (8.4) 



P 
j 



w + w 

2vr_ 

£ ~ £ 
2m ~ 
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We substitute these expressions back into (8.2) and obtain 

S[u] = J d 2 x [ e (u,) -£(-«>)] + J W ^#)J. (8.5) 

We assume now that the integration variable of the path integral is the displace- 
ment field u. The fields w, w in (8.5) are not independent but are related to u 
by (8.4). Deriving (8.5) we neglected the fluctuations around the saddle point 
as well as the changes in the measure of path integration. The corrections due 
to the neglected terms will be of the higher order in field gradients and are not 
essential for our calculation of the leading term of the EFP. 

One can now use (8.5) in the path integral formulation of the quantum hy- 
drodynamics (2.1) where the integration is taken over all configurations of the 
displacement field u. We note here that in the Galilean invariant systems the 
"kinetic term" of the hydrodynamic action (3.1) is fixed. The Lagrangian of 
the action (8.5) is a more complicated function of dtu. In this case we derived 
it using the fact that fermions are free. In a more general problem of interacting 
particles the derivation of the hydrodynamic action requires the solution of the 
dynamic many body problem. 

The Lagrangian of (8.5) is given by 

1 r w rlk 

L = -^ w ~^ t £ H " + / _ 2^ £{k) - (8 - 6) 

Now, we calculate the EFP for lattice fermions using the the Lagrangian (8.6) 
with the Eq. (8.4), and the results of Appendix A. First, we calculate 

dL= £ -±l dp+ ^ dj (8.7) 

and 

Lj = v = ^, (8.8) 

L P = ( 8 .9) 

^ = ^ (B-10) 
L - 271 £ '~"' (8 11) 

where e' = 1 1| | k=w and e' = ||| | fc=j _ . In terms of the density and velocity 

w = irp + iv. (8.13) 
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We emphasize that because of the absence of the Galilean invariance (the dis- 
persion is not k 2 /2m), the current is not pv but can be found from (8.4,8.8). 
We obtain then an interesting relation 



K= y /L pp L jj -L 2 pj = ir. (8.14) 

This relation can be considered as a special property of free fermions (with 
arbitrary dispersion!). Namely, one can trace the origin of (8.14) to the fact 
that each fermion occupies a fixed volume in the phase space. 7 

Let us now apply Riemann's trick to equations (8.3). Namely, we inter- 
change independent and dependent variables so that x = x(w,w) and t = 
t(w, w). Then using (8.14, 8.8-8.12) we re-write (8.3) as 

d w (x-ie't)=0 (8.15) 

and its complex conjugate. The most general solution of (8.15) is given im- 
plicitly by 

w = F(z), (8.16) 
where F(z) is an arbitrary analytic function and 

z = x + ie'{w)t. (8.17) 

The equations (8.15, 8.17) reduce to equations for free continuous fermions if 

e(k) = k 2 /2. 

Similarly to free continuous fermions, the EFP boundary conditions specify 
the form of the function F(z) in (8.16). We define 

w — irp „ „ 

- = wrr <8 - 18) 

where we assume lattice spacing to be 1 so that the maximal density of fermions 
on the lattice is 1. Then the EFP solution is given by F(w) such that 

z 

s'mw = sin 6 (8.20) 

Vz 2 - R 2 



We immediately find at large distances (z — > zq = x + ie' t — ► oo) 

R 2 



w — npo 2(1 — p) — jtan^ (8.21) 
2z n 



7 For Calogero-Sutherland model (7.4) we have k = nX. It means that the volume of the phase space per 
particle is changed by the factor of A which reflects the fractional statistics of particles. 



d P0 S = 2tt(1 - p)R 2 tan ^ u ^ , (8.23) 
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and using (4.9) and (8.13) we extract 

a = — (1 -p)i? 2 tan 0. (8.22) 

Z7T 

We substitute this expression in (A.7) and obtain 

7t(P0 - p) 

2(1 -p) 
so that 

S DFP = -4(1 - pfB? In cos . (8.24) 

In particular for p = we obtain for the Emptiness Formation Probability 

5 BJ?P = -4i? 2 lncos^ = -4^111 cos ^. (8.25) 

This result is the exact first term in l/R expansion (B.10). The next term is 
\ In R. It is interesting to note that as it should be, in the limit po — p <C 
po the result (8.24) reproduces the bosonization result (4.10). Also, in the 
limit po,p <C 1 the results (8.24,8.25) reproduce the corresponding results 
(6.17,6.18) for free continuous fermions. 



9. Conclusion 

We showed on the example of the Emptiness Formation Probability that one 
can obtain some exact results for correlation functions using a collective hy- 
drodynamic description. Moreover, the nonlinear hydrodynamic description 
might work and produce non-perturbative results even in cases where the lin- 
earized hydrodynamics fails. We obtained the leading term of the asymptotics 
of the EFP using the instanton approach. Although the EFP is a property of 
the ground state of a quantum many body system, the hydrodynamic approach 
can also be used to study dynamics of quantum systems. It is also not limited 
to integrable systems. 

We considered few simple systems which are in the Luttinger liquid phase 
at zero temperature. The effects of finite temperature and finite gap in the spec- 
trum of excitations can also be considered using the hydrodynamic approach 
in the limit when temperature is very low and the gap is very small [10, 12]. In 
these limits one obtains a crossover between the Gaussian decay of the EFP at 
intermediate R to the exponential decay at very large R. 

There are also many questions which are left open. First, it would be nice to 
obtain the results for the EFP in other integrable systems such as bosons with 
delta-repulsion and XXZ spin chains. Some results for the EFP have already 
been obtained using other methods (see Appendix B), while, e.g., the EFP for 
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a XXZ spin chain in the presence of magnetic field is still not known. The 
application of the hydrodynamic approach developed here to these systems is 
straightforward and reduces the problem to the problem of finding asymptotics 
of the solution of a system of classical equations with proper boundary condi- 
tions. However, these classical equations are complicated and we haven't yet 
obtained analytical results for their asymptotics. 

The second open question involves the corrections to the leading hydrody- 
namic approximation that we used in this lecture. There are two sources of 
such corrections: the gradient corrections to the "classical" hydrodynamic ac- 
tion and quantum fluctuations around the classical saddle point. E.g., if one is 
interested in the next to leading terms of the asymptotic expansion of the EFP 
in l/R one needs to include these gradient corrections (which make hydro- 
dynamics dispersive) and quantum fluctuations. Especially interesting would 
be to obtain the power law pre-factor of the EFP (or In R term in the asymp- 
totic expansion). The hydrodynamic calculation might shed some light on a 
possible universality of the exponent of this pre-factor. 

Needless to say, that although we focused in this lecture on the particu- 
lar correlation function (the EFP), the use of hydrodynamic approach is much 
broader. In particular, the nonlinear hydrodynamics is important to capture a 
lot of important nonlinear phenomena which disappear in the linear approxi- 
mation. 
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Appendix: Hydrodynamic approach to non-Galilean invariant 

systems 

The hydrodynamic description of Sec. 1.3 and the approach of Sec. 1.5 can be easily gener- 
alized to non-Galilean invariant systems. We briefly list here relevant formulas leaving the de- 
tails of (straightforward) calculations to the interested reader. We assume that a one-dimensional 
compressible liquid can be described by the partition function (2.1) where the functional integra- 
tion is taken over all configurations of the displacement field u which defines the hydrodynamic 
density and current by (3.3) and where the Euclidian action (2.2) can be written as 



The Lagrangian density L(p,j) is a function of p and j. As in the main text we neglect the 
gradient corrections to the action (A.l). In the case of non-Galilean invariant systems these 




(A.1) 
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gradient corrections will depend on gradients of current in addition to gradients of density. 
The case of Galilean invariant systems considered in the main text (3.1) can be obtained from 
the general one assuming the particular Galilean invariant form of the Lagrangian L(p,j) — 
f/2p + pe(p). 

Variation of (A.l) with respect to u gives a generalized Euler equation which together with 
the continuity equation reads 

d t p + d x j = 0, (A.2) 
d t Lj - d x L p = 0, (A.3) 

where Lj means dL/dj etc. 

If the deviation from the ground state is small, we can expand the action (A.l) around p = 
po, j — 0. We assume that L*v = 0, where superscript (0) means that the derivative is 
calculated at the equilibrium values p = po and j = 0. We obtain for the action in harmonic 
approximation 

S = J (fx [Lf)u 2 t + 4°>i£] . (A.4) 
We introduce a complex space-time coordinate 

2o = x + iv s0 t, (A. 5) 

where 



V S = 



(0) 

(A.6) 



is the sound velocity at equilibrium. Let us define the coefficient a through the asymptotics of 
an EFP solution at zo — > oo by (4.7). The asymptotics of the current and density are given by 
(4.9). 

Our main formula (5.4) for the variation of the EFP instanton action with respect to the 
background density becomes 

d P0 S — 2ttk (a + a), (A.7) 

where kq is given by 

«o = y/L?}L?>. (A.8) 

The coefficient ko is related to the compactification radius i? C om P of bosons in the bosonization 
procedure k = (2nR comp ) 2 . For free fermions n = n which corresponds to R comp = 

l/v 7 !^. 



Appendix: Exact results for EFP in some integrable models 

For the sake of reader's convenience in this appendix we list some results obtained for the 
Emptiness Formation Probability P(R) in integrable one-dimensional systems. We present the 
results for S = — In P(R) which should be compared with the instanton action S op t used in 
the main text. 



Free continuous fermions Let us denote 

S EE TVp R. (B.l) 

We use the fact that the ground state wave function (more precisely |*| 2 ) of free fermions 
coincides with the joint eigenvalue distribution of unitary random ensemble. For the latter the 
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probability of having no eigenvalues in the range 2R of the spectrum was obtained in [8] (see 
also [9]). First few terms of the expansion in 1/R are 

S=±s 2 + \]ns- ^In2 + 3C'(-1)) +0( S - 2 ). (B.2) 

Calogero-Sutherland model The Calogero-Sutherland model [14, 15] (rational ver- 
sion, known also as Calogero model) with iV-particles is defined as 8 

n, 1^ 2 1 A(A-l) 

7 = 1 i,k=l:ijtk 



where pj = —id/dxj is the momentum operator of j-fh particle and A is a dimensionless 
coupling constant. The wave function of the ground state is proportional to Ylj < k( x i i — x k) • At 
A = 1 we have free fermions, while in the case of general A the model (B.3) describes particles 
with fractional statistics. Using the form of the ground state wave function and thermodynamic 
arguments [9] one obtains 

S= ^s 2 + (l-A)s + 0(lns). (B.4) 

or defining 

s = v^XirpoR. (B.5) 

and 

A l/2 _ A -l/2 

qo = (B.6) 

we have 

S= is 2 -2a ( ,s + 0(lns). (B.7) 

The notation o?o originates from the conformal field theory with central charge c = 1 — 24a 2 , 
which is known to be related to the Calogero-Sutherland model [16]. As far as I know the 
coefficient in front of In s term of the expansion is not known for the general A. However, 
A = 1/2, 1, 2 correspond to random matrix ensembles where the full asymptotic expansion is 
known (see below). In those cases the coefficient of Ins is 1/8, 1/4, 1/8 respectively. The 
natural guess is that 

S = is 2 - 2a s "o) Ins + 0(1). (B.8) 

Random matrices For Random Matrix ensembles with f3 = 1, 2, 4 the joint eigenvalue 
distribution is proportional to Yii<j l z * ~ z i \^ ■ The full asymptotic expansion of the quantity 
Ep(0, 2R) corresponding to the EFP P(R) was obtained using properties of Toeplitz determi- 
nants [8, 9]. The first few terms of these expansions are given by 

Sx=i/2 = i(^oi?) 2 + i(^ofi) + iln(vTp ()J R)-^ln2-|c'(-l) + 0(l/s), 

Sx=i = i(7rpoi?) 2 + iln(vTp i?)-^ln2-3C'(-l) + 0(l/ S 2 ), (B.9) 

Sa= 2 = (vTp OJ R) 2 -(vTpo J R) + iln(7Tp„i?) + ^ln2-^C'(-l) + 0(l/s). 

o o Z 



8 To prevent particles running to infinity we either add a harmonic potential to (B.3) or put particles on a 
circle of a large radius. 
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Here we used A = (3/2 = 1/2, 1, 2 instead of (3. Using notations (B.5,B.6) we can summarize 
the first three terms of (B.9) in a compact form (B.8). 

Free fermions 011 the lattice For non-interacting one-dimensional fermions on the lat- 
tice (and the corresponding XY spin chain) the asymptotic behavior of EFP was derived in [17] 
using the Widom's theorem on the asymptotic behavior of Toeplitz determinants. Introducing 
the Fermi momentum &f = npa and using units in which the lattice spacing is 1 we have 



S = -4R 2 In cos — + - In 
2 4 



2i?sin^ 
2 



L2 ln2 + 3C'(-1)) +0(R- 2 ). (B.IO) 



In the continuous limit &f — > the (B.IO) goes to its continuous version (B.2). 

Bosons with delta repulsion The model of bosons with short range repulsion is de- 
scribed by 

1 N 

n =2^ti +g X) 8(xj-x k ), (B.ll) 

3=1 l<j<k<N 

where g is a coupling constant. It is integrable by Bethe Ansatz [18]. It was derived (conjec- 
tured) in Ref. [19] that the leading term of the EFP is 

S=\{KRf[l + I{g/K)\, (B.12) 
where K is the Fermi momentum in the Lieb-Liniger solution [18] and 

m = J_ f 1 y d y f 1 *Ji log ( x 2 + (y + z) 2 \ (B 13) 

The limit I(x — > oo) = corresponds to the free fermion result (B.2) (Tonks-Girardeau gas), 
while the limit I(x — » 0) = 1 is the result for the EFP of free bosons. 

XXZ model The Hamiltonian of an XXZ model is given by 

H = J ^2 [0k<Tfc+i + o- y k a y k+1 + Aa z k a z k+1 ] , 

k 

where the sum is taken over the sites of a one-dimensional lattice and a x ' y ' z are Pauli matrices. 
Let us parametrize the anisotropy as A = cos irv. Then for the EFP we have [20, 21] 

P(n) -An-'C-" 2 , (B.14) 

as n = 2R — > oo, where 

r _ T 2 (l/4) f rdt sinh 2 ^^ \ 

ttV^ P l Jo t cosh(2^)sinh(t)/ {tS - n) 

and the exponent 7 was conjectured in [21] to be 

1 2 

7= h —A r- (B.16) 

' 12 3(l-i/) 
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